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Stellar companions to exoplanet host stars: 
Lucky Imaging of transiting planet hosts* 
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ABSTRACT 

Observed properties of stars and planets in binary/multiple star systems provide clues 
to planet formation and evolution. We extended our survey for visual stellar compan- 
ions to the hosts of transiting exoplanets by 21 stars, using the Lucky Imaging tech- 
nique with the two AstraLux instruments: AstraLux Norte at the Calar Alto 2.2-m 
telescope, and AstraLux Sur at the ESO 3.5-m New Technology Telescope at La Silla. 
Typically a sensitivity to companions of magnitude difference Az' ~ 4 is achieved at 
angular separation p = 0.5" and Az' > 6 for p = 1". 

We present observations of two previously unknown binary candidate companions, 
to the transiting planet host stars HAT-P-8 and WASP-12, and derive photometric 
and astrometric properties of the companion candidates. The common proper mo- 
tions of the previously discovered companion candidates with the exoplanet host stars 
TrES-4 and WASP-2 are confirmed from follow-up observations. A Bayesian statistical 
analysis of 31 transiting exoplanet host stars observed with AstraLux suggests that 
the companion star fraction of planet hosts is not significantly different from that of 
solar-type field stars, but that the binary separation is on average larger for planet 
host stars. 
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1 INTRODUCTION 

About half of solar-type stars in our neighbourhood are part 
of a binary or multip le system (|Duquennov fc Mavor|[l99ll ; 
iRagh avan e t al. 2010). Understanding how a secondary star 
affects the formation and evolution of planets in the system 
is therefore of high importance for an estimate of the overall 
occurrence of planets in our Galaxy. 

A close stellar companion is expected to affect 
planet formation in several ways, e.g., by heating 
and truncating the ci rcum st ellar protoplanetary disc 
llArtv mowicz fc Lubowl 1 1994 lArmitage. Clarke, fc ToutJ 
1 19991 : iNelsonl |2000D . or by increasing the relative veloc- 



* Based on observations collected at the 2.2 m telescope at Calar 
Alto Observatory, and on observations made with the 3.5 m ESO 
NTT at La Silla Observatory under programme ID's 082.C-0084 
and 083.C-0145. 
f E-mail: bergfors@mpia.dc 



ities of the planete simals (|Heppenheime 3 11974 Il978l ; 
IWhitmire et al.l 19981) . A sec ondary star may either stim- 
ulate (|Bossl l2006h or hinder (jNelsonl l200d : iKlev fc Nelson] 
|200ST ) planet formation. The occurrence and properties of 
planets formed in binary systems may also provide a way to 
discriminate between the two most wid ely supported planet 
formation models: core accretion (e .g.. | Pollack et al.lll996f) 



and g ravitational instability (e.g., iBossf 1997 : Mayer et al.l 
2002). For instance, formation by gravitational instability 
(GI) requires a massive disc for fragmentation to occur 
(M > O.IMq). iMaver et all (|2005h found that planet for- 
mation by GI in massive discs is not significantly different 
from formation around a single star if the binary separation 
is > 120 AU, while at close binary separations (< 60 AU) 
the high temperatures caused by shock heating effectively 
suppress fragmentation. Formation by core accretion, on 
the other hand, does not require a massive disc, and planets 
may form closer to their parent star without being much 
affected by a close binary companion. Therefore, a trend 
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of giant planet frequency with binary star separation, i.e., 
fewer giant planets in binary systems closer than 100 AU 
than in wider binaries or single s ystems, might poin t to GI 
as a main formation mechanism (|Maver et alj|200a ). 

Observations show that the frequency by which giant 
planets are formed in close binary stellar systems appears to 
be slightly lo wer than around single stars. lEggenberger et al.l 
(|200a . l201lh found that giant planets are more common m 
single-star systems than in binaries separated by 35-100 AU. 
Nevertheless, surveys of multiplicity among planet host stars 
show that many exoplanet host stars are part of a binary or 
multiple system. We now know of ~ 50 binary or multiple 
systems where the planet (s) belong to one of the stars (S- 
type orbit), and Kepler observations recently revealed the 
first known planets that tran sit both stars in a close binary 
syste m, i.e., in a P-type orbit (iDovle et al .l201ll ; IWelsh' et ajj 
l2012h . System characteristics such as binary separation to- 
gether with properties of the planets (orbital period, mass, 
eccentricity, etc.) and any differences compared to the prop- 
erties of the single-star planetary systems provide important 
constraints on planet formation and dynamical system evo- 
lution. Among suggested correlations between stellar and 
planetary properties in binaries, two appear significant: The 
most massive planets in short-period orbits belong to stars 
in binary sy stems, and so do the planets with the highest ec- 
centricities (iDesidera fc Barbieril [20071 ; ISozzetti fc Desideral 
l2010l ; lTamuz et al.ll200Sl )~ 



A significant fraction of exoplanets to date have been 
discovered using radial velocity measurements, and such 
systems have also b een the subjec t of several multiplic- 



20061; Eggenberger et all 2007 


; Mugrauer et all 120071: 


Roberts et al.ll201ll; iMason et al. 


l201ll; iGinski et al.|[2012). 



However, the discovery rate of transiting exoplanets 
(TEPs) has increased tremendously over the last couple 
of years, mainly due to the success of ground-based tran - 
sit searches s uch as SuperWAS P (|Pollacco et al.l 120061 ) 
and HATNet (|Bakos et al.l |2004). as well as space-based 
programs such as C 0R0T (|Baglin et al.l l2006h and Kepler 
l|Borucki et al.ll201Ch . Transiting exoplanets are unique in 
the way that properties such as mass and radius can be 
measured, and a variety of additional physical parameters 
such as true mass, mean density and surface gravity can 
thus be derived, potentiall y providing information on 
planet formation (see, e.g.. iMazeh. Zucken fc Pont 20051 ; 
iTorres. Winn, fc Holmanl |200S| ; ISouthwortbl |2009| . (2010t ). 
An unresolved faint secondary star within the photometric 
aperture, whether it is a bound companion or chance 
alignment, contributes a constant flux offset to the transit 
light curve and affects the accuracy with which stellar 
and planetary parameters can be derived. High resolution 
imaging of TEP hosts have shown that the presence of a 
(projected) nearby star may require a correction of derived 
stellar and p lanetary parameters between a few to seve ral 
tens percent l|Daemgen et al.ll2009l ; iBuchhave et alj|201ll ). 

Another important aspect of transiting planets is that 
they allow for a measurement of the projected alignment 
between the orbital plane of the planet and the rotational 
plane of the star through the Rossiter-McLaughlin effect 
(e.g. lOueloz et al.l l2000t ). Studies using this effect have 
shown that many transiting hot Jupiters are significantly 
misaligned with respect to the stellar spin (e.g. I Winn et al.l 



l201dt) . which implies that three-body mechanisms such as 
Kozai migration (e.g. iFabrvckv fc Trem ainc 2007]) may be 
responsible for the formation of those systems, rather than 
classical orbital migration (e.g. iLin et al1ll996l h Since the 
Kozai mechanism requires the presence of a wide stellar com- 
panion, it follows that searching for binarity among TEP 
hosts could provide important clues for the formation of hot 
Jupiter systems. 

In this paper we present high resolution Lucky Imag- 
ing observations of 21 TEP host stars using the two As- 
traLux instruments at the 2.2 m telescope at Calar Alto and 
at NTT at La Silla. The observ ations complement the sam- 
ple of TEP hosts presented by iDaemgen et all ([2009) , and 
are described in Sect. 2 together with the methods used 
for obtaining relative astrometry and photometry. All the 
transiting exoplanets in our survey transit only one star, 
in a short-period orbit. In Section 3 we derive magnitudes, 
(i — z) colours, and photometric spectral types and distances 
to the companion candidates from the photometric observa- 
tions in SDSS i! - and z'-band and known spectral types of 
the planet host stars. Evidence of physical companionship 
of the planet host stars and the companion candidates is 
investigated, and the results for each individual target are 
compared with previously published astrometric and pho- 
tometric data if available. We estimate the probability of 
chance alignment and perform a Bayesian analysis of our 
complete observed sample of TEP hosts with AstraLux. The 
results are summarized and discussed in Sect. 4. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations with AstraLux 

The 21 transiting exoplanet host stars were observed within 
four different observing runs with the two AstraLux Lucky 
Imaging instruments. Most of the TEP hosts in the sur- 
vey are located in the northern sky and were therefore ob- 
served with AstraLux Norte at the 2.2 m telescope at Calar 
Alto observatory in October-November 2009, with follow-up 
astrometric observations in November 2011. The southern- 
sky targets were observed in November 2008 and April 2009 
with the AstraLux Sur visitor instrument mounted to the 
ESO 3.5 m New Technology Telescope (NTT) at La Silla 
(see iHormuth et~al"1 120081 ; [Hippler et at] 120091 . for details 
on the AstraLux instruments). Lucky Imaging is a way to 
limit the effects of atmospheric turbulence by taking a large 
number of very short integrations from which only the least 
distorted few percent of the frames are selected. These are 
shifted and added to produce the final image, yielding al- 
most diffraction-limited resolution. The full AstraLux field 
of view (FoV) in the final resampled frames is ~ 15.7" x 15.7" 
for AstraLux Sur, and « 24" x 24" for AstraLux Norte. The 
individual exposure time was either 15 ms or 30 ms, depend- 
ing on the target brightness and observing conditions. The 
shorter integrations were achieved by decreasing the FoV 
and reading out only a subframe of the detector. In order 
to match a total integration time of 300 s, the number of 
integrations was set to 20 000 or 10 000 respectively. 

Each target was observed in SDSS i'- and z'-filter. As- 
trometric reference stars in the open cluster NGC 3603, the 
Orion Trapezium and/or the globular cluster 47 Tuc were 
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Table 1. Transiting oxoplanet host stars with no observed com- 
panions. 



Planet Host Instrument 



FoV Obs. Date (UT) 



2 3 4 5 

Angular Separation ["] 

Figure 1. Average sensitivity of the AstraLux Norte observa- 
tions in October-November 2009. The squares mark the compan- 
ion candidates (see Table [2j- The dashed line represents the typ- 
ical 5<r detection limit. 



observed each night for calibration of the detector rota- 
tion and pixel scale. Field rotation and plate scale were 
derived using custom IDL routines to compute separation 
and position angle of the astrometric reference stars pair- 
wise in several frames and compare to values derived from 
e.g. HST/WFPC2 archive images. The plate scale for the 
observations with AstraLux Norte was 23.43 ± 0.06mas/px 
with the detector rotated 0.06 ± 0.02° east of north in 
October-November 2009 and plate scale 23.74 ±0.05 mas/px 
in November 2011 with a detector rotation of 1.6 ± 0.2° 
to the east. For the AstraLux Sur instrument we derived 
a plate scale of 15.373 ± 0.002 mas/px and detector rota- 
tion 1.7 ± 0.3° to th e west of north in November 2008 (see 
iBergfors et ai1l2010h . and 15.245 ± 0.006 mas/px, 1.4 ± 0.2° 
to the west in April 2009. 



2.2 Photometry and astrometry 

In the sample of 21 TEPs, we found candidate companion 
stars from visual inspection of the reduced Lucky Imaging 
frames to 7 stars: t he previously known compan ion can- 
didates to WASP- 2 (ICollier Cameron et all [20071). TrES-2 
and T rES-4 l|Daemgen et al.ll2009T ). HAT-P-7~ (|Narita et all 
2010), and new companion candidates to HAT-P-8, WASP- 
12 and XO-3. Stars for which no companions were detected 
within the FoV of AstraLux in these observation runs are 
listed in Table [1] The stars for which we observed a candi- 
date companion are listed in Table [2] 

Most of the companion candidates reside close to the 
primary stars, within the PSF-wings (HAT-P-7 and XO-3 
being the exceptions). For these stars we performed relative 
photometry and astrometry of the companion c andidates 
using mainly the IRAF allstar |Todvlll986l , Il995t ) task for 
PSF-fitting. The PSF was built from the primary star in 
an image where the secondary star had been removed (see 
iDaemgen et al.ll2009l . for more details on the procedure). 
Aperture photometry with IRAF apphot was used to de- 
termine the properties of the wide companion candidates to 
XO-3 and HAT-P-7. The astrometric and photometric prop- 
erties in Tables 131 II are averaged measurements of the final 



HAT-P-2 


AstraLux 


Norte 


12.00 X 12.00 


1 


Nov 


2009 


HAT-P-11 


AstraLux 


Norte 


12.00 X 12.00 


29 


Oct 


2009 


HAT-P-13 


AstraLux 


Norte 


12.00 x 12.00 


2 


Nov 


2009 


HD 149026 


AstraLux 


Norte 


12.00 x 12.00 


1 


Nov 


2009 


HD 209458 


AstraLux 


Sur 


7.87 x 7.87 


12 


Nov 


2008 


HD 80606 1 


AstraLux 


Norte 


12.00 x 12.00 


30 


Oct 


2009 


WASP-3 


AstraLux 


Norte 


12.00 x 12.00 


30 


Oct 


2009 


WASP-4 


AstraLux 


Sur 


15.74 x 15.74 


10 


Nov 


2008 


WASP-5 


AstraLux 


Sur 


15.74 x 15.74 


10 


Nov 


2008 


WASP-7 


AstraLux 


Sur 


7.87 x 7.87 


10 


Nov 


2008 


WASP-13 


AstraLux 


Norte 


12.00 x 12.00 


30 


Oct 


2009 


WASP-15 


AstraLux 


Sur 


7.83 x 7.83 


12 


Apr 


2009 


XO-4 


AstraLux 


Norte 


12.00 x 12.00 


1 


Nov 


2009 


XO-5 


AstraLux 


Norte 


12.00 x 12.00 


31 


Oct 


2009 



M Visual binary star (e.g., iNaef et alj|200lh . outside of 
AstraLux field of view. 

images using the Lucky Imaging combination of the best 
5% and 10% of the individual integrations. For follow-up 
astrometry in November 2011, only z'-band measurements 
were used. The separation and position angle error bars are 
propagated from la uncertainties in the positional measure- 
ments, and take into account systematic errors (i.e., plate 
scale and detector orientation, see Sect. 2.1). The uncer- 
tainty in detector rotation is usually the dominant error 
source in the relative astrometry. The error bars in magni- 
tude differences are propagated from the photometric errors 
estimated with IRAF. 

Figure [T] shows the typical 5cr detection limit for the ob- 
servations with AstraLux Norte in October-November 2009. 
At separation p — 0.5" we can typically detect a companion 
4 magnitudes fainter than the primary in Az' , and 6 mag- 
nitudes fainter at p = 1". A similar plot repr esenting obser- 
vation s with AstraLux Sur can be found in IBergfors et al.l 
|20ld) . The Amag as a function of angular separation was 
estimated by pairwise subtraction of a set of images of single 
stars with similar peak flux and FWHM from each other and 
measuring the noise level in concentric circles at increasing 
separations from the centre. 



3 RESULTS 

3.1 Properties of the stellar companion candidates 

The spectral types of the companion candidates are derived 
using primary star spectral types from literature and assum- 
ing that the additional flux from the much fainter secondary 
sta rs do not affect these spectra l types. We use stellar SEDs 
by iKraus fc Hillenbrand! (|2007h for the known TEP host 
spectral type in combination with 2MASS photometry in 
the JHK bands (|Cutri et al.| [2003) to find the average dis- 
tance modulus and SDSS i- and z-band magnitudes. The 
measured magnitude difference between the primary and 
secondary stars in i' and z then provides apparent mag- 
nitudes and (i — z) colour for the companion candidate. 
While we do not explicitly correct for the transformation 



© 2012 RAS, MNRAS 000, [TH] 



4 C. Bergfors et al. 




Figure 2. AstraLux Norte z 1 — band observations of the candidate binaries (a) HAT-P-8 and (b) WASP-12 in October 2009. The images 
are shown in a square root scale with north up and east to the left. 



Table 2. Relative photometry and astrometry for the companion candidates. 



ID 


Date of obs. 


Instrument 


FoV 


Separation 


Position Angle 


Az' 


Ai' 








["] 


["] 


[°] 


[mag] 


[mag] 


HAT-P-7 1 


30 Oct 2009 


AstraLux Norte 


12.00x12.00 


3.82±0.01 


90.4±0.1 


6.89±0.07 


7.94±0.21 




9 Nov 2011 


AstraLux Norte 


12.00x12.00 


3.86±0.07 


89.9±0.8 


6.87±0.07 


7.58±0.07 


HAT-P-8 


29 Oct 2009 


AstraLux Norte 


12.00x12.00 


1.027±0.011 


137.3±0.4 


6.68±0.07 


7.34±0.10 


TrES-2 2 


29 Oct 2009 


AstraLux Norte 


12.00x12.00 


1.085±0.006 


136.1±0.2 


3.48±0.06 


3.73±0.03 


TrES-4 2 


30 Oct 2009 


AstraLux Norte 


12.00x12.00 


1.550±0.007 


359.9±0.2 


4.19±0.05 


4.57±0.05 


WASP-2 3 


13 Apr 2009 


AstraLux Sur 


15.66x15.66 


0.761±0.009 


103.5±0.2 


3.62±0.05 






29 Oct 2009 


AstraLux Norte 


12.00x12.00 


0.739±0.024 


104.0±1.3 


3.64±0.04 


4.17±0.03 




9 Nov 2011 


AstraLux Norte 


12.00x12.00 


0.744±0.013 


104.6±0.7 






WASP-12 


30 Oct 2009 


AstraLux Norte 


12.00x12.00 


1.047±0.021 


249.7±0.8 


3.79±0.10 


4.03±0.07 




8,9 Nov 2011 


AstraLux Norte 


12.00x12.00 


1.043±0.014 


249.9±0.5 






XO-3 


30 Oct 2009 


AstraLux Norte 


12.00x12.01 


6.059±0.047 


296.7±0.3 


8.22±0.23 


8.57±0.24 



W The October 2009 observations were fi rst published in |Nj3xita_et_al 1 l|201Cf l M Companion was disco vered bv | D aemgen et al . (2009) 
l 3 l Companion was discovered by [Collier Cameron et al. I ll2007l l and has also been observed by iDaemgen et alTl200*9fi 



between SDSS i',z' to i, z, the difference, when calculated 
from the photometric transformation equations (SDSS web- 
page), is small. The magnitude difference between the two 
photometric systems is less than (i — i 1 ) = 0.05 mag for stars 
bluer than (r — i) « 1.5, and even less in (z — z'). Accord- 
ing to the SEDs of lKraus fc Hillenbrand! (|2007f) . this (r - i) 
colour corresponds to spectral types M4 or earlier. None of 
our derived spectral types are later than M4, and we conser- 
vatively assume error bars of ±0.1 for the derived i, z mag- 
nitudes of the primary and secondary stars so as to include 
the transformation between photometric systems as well as 
uncertainty in the photometric measurements (±0.2 for the 
very faint companion candidate to XO-3). The secondary 
star spectral types and primary and secondary photometric 
distances w ere estimated from the deriv ed (i — z) colour and 
the SEDs o f K raus fc Hillenbrand! (|2007t ). Interstellar extinc- 



tion and stellar metallicity were not considered in these es- 
timates but may affect the photometric distances, as well 
as the colours and spectral types. The photometric distance 
estimates are therefore only indicative and not hard limits. 
The components' spectral types, apparent magnitudes, pho- 
tometric distances and the secondary stars' (i — z) colours 
are listed in Table [3] 



3.2 Notes on individual systems 

Companion candidates to 7 transiting exoplanet hosts were 
found from these observations. Of these, the three TEP 
hosts WASP-2, TrES-2 and TrES-4 and their candidate 
companions had been observed with AstraLux previously 
l|Daemgen et all 120091 ). a nd the cand i date companion to 
HAT-P-7 is discussed in iNarita et all (^OloT ). The widely 
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Table 3. Photometric properties of the primary stars and the companion candidates. 



Planet Host 


SpT(A) 




SpT(B) 


m z (A) 


m ; (A) 


m z (B) 


m ; (B) 


(i - z) B 


d(A) 




d(B) 










[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


be] 




[PC] 


HAT-P-8 1 


F8V 


Ml V ... 


M3.5V 


10.20±0.1 


10.18±0.1 


16.88±0.1 


17.52±0.1 


0.64±0.14 


150±20 


710 . 


.. 350 




F5V 


Ml V ... 


M3.5V 


10.08±0.1 


10.03±0.1 


16.76±0.1 


17.37±0.1 


0.61±0.14 


190±20 


670 . 


.. 330 


TrES-2 


GOV 


KOV 


... MOV 


ll.lliO.l 


ll.llztO.l 


14.59±0.1 


14.84±0.1 


0.25±0.14 


220±20 


670 . 


.. 270 


TrES-4 


F8V 


K4.5 V ... 


M1.5V 


11.57±0.1 


11.55±0.1 


15.76±0.1 


16.12±0.1 


0.36±0.14 


290±30 


790 . 


.. 380 


WASP-2 2 


KOV 


Ml V ... 


M3.5V 


H.lOiO.l 


11.20±0.1 


14.74±0.1 


15.37±0.1 


0.63±0.14 


140±10 


260 . 


.. 130 




K2V 


M1.5V 


... M4V 


11.19±0.1 


11.33±0.1 


14.83±0.1 


15.50±0.1 


0.67±0.14 


120±10 


240 . 


.. 110 


WASP-12 


GOV 


KOV 


... MOV 


11.41±0.1 


11.41±0.1 


15.20±0.1 


15.44±0.1 


0.24±0.14 


250±30 


890 . 


.. 350 


XO-3 3 


F5V 


GOV ... 


M2.5V 


9.91±0.1 


9.86±0.1 


18.13±0.2 


18.43±0.2 


0.30±0.30 


170±20 


5470 


...860 



M A primary star spectral type of F8V is suggested by|joncs & Sleep (201(3), while we find F5 V an equally good fit to the SEDs of 

iKraus fc Hillenbran 3 J2007I) . We therefore list both alternatives here. I 2 ' The primary star is of spectral type Kl V. 
IKraus fc Hillenbrand! ^POTt i provides spectral type models only for KOV and K2V, and we choose to present these two alternatives 
rather than interpolate between spectral types. PI Here we assume that the companion candidate is a main sequence star when deriving 

spectral type and distance, however see Sect. 3.2.6 



Table 4. Summary of astrometric measurements. 



Planet Host Date of Obs. Separation Pos. Ang. Ref. 

H [°] 



HAT-P-7 


6 Aug 2009 


3.88±0.01 


89.8±0.3 


[1] 




30 Oct 2009 


3.82±0.01 


90.4±0.1 


[1],[2] 




9 Nov 2011 


3.86±0.07 


89.9±0.8 


[2] 


HAT-P-8 


29 Oct 2009 


1.027±0.011 


137.3±0.4 


[2] 


TrES-2 


May 2007 


1.089±0.008 


135.5±0.1 


[3] 




29 Oct 2009 


1.085±0.006 


136.1±0.2 


[2] 


TrES-4 


Jun 2008 


1.555±0.005 


359.8±0.1 


m 




30 Oct 2009 


1.550±0.007 


359.9±0.2 


[2] 


WASP-2 


Nov 2007 


0.757±0.001 


104.7±0.3 


[3] 




13 Apr 2009 


0.761±0.009 


103.5±0.2 


[2] 




29 Oct 2009 


0.739±0.024 


104.0±1.3 


[2] 




9 Nov 2011 


0.744±0.013 


104.6±0.7 


[2] 


WASP-12 


30 Oct 2009 


1.047±0.021 


249.7±0.8 


[2] 




8,9 Nov 2011 


1.043±0.014 


249.9±0.5 


[2] 



W iNarita et all j201Ct) ' 2 1 This paper I 3 ) iDaemgen et al.l d200gf) 



separated candidate companion to XO-3 is likely to be a 
physically unrelated background object if it is a main se- 
quence star, although the possibility of a coeval white dwarf 
companion can not be ruled out from these observations. 
The two faint objects at separations of p ~ 1" to HAT-P-8 
and WASP-12 are previously unknown, plausibly bound stel- 
lar companions. Figure[2]shows the AstraLux Norte z'— band 
observations of these two systems obtained in October 2009. 
All astrometric measurements including previous observa- 
tions are summarised in Table [4] 



3.2.1 HAT-P-8 

The transitin g exop lanet HAT-P-8 b was discovered by 
Latham et al l (|2009h . It is a slightly inflated planet with 
mass M p = 1 52t^ gMj and radius R p = 1.50±£ ^Rj 
|Latham et al.l 120091 '). The host star spectral type is only 
giv en as F in the discove ry paper, but is referred to as F8 V 
by IJones fc Slee^l (|2010f i. We find from the 2MASS JHK 
photometry and colours that a spectral type of F5 V also 



fits the SEDs of lKraus fc Hillenbrand! (|2007l ), and list both 
alternatives in Table until a more precise spectral classi- 
fication can be made. Although the assumed primary spec- 
tral type affects the derived primary i',z' magnitudes and 
thereby the colour and spectral type of the companion candi- 
date (see Sect. 3.1,) we find for both alternatives of primary 
spectral type that the stellar candidate companion is likely 
to be of spectral type M2 to M4 from the (i — z)-colours. 



3.2.2 TrES-2 

The companio n can didate was first discovered by 
IDaemgen et~aH (|2009l ) from AstraLux Norte observa- 
tions in May 2007. The observations presented here 
took place in October 2009. We measure a separation 
of ARA = 0.752" ± 0.007", ADec = 0. 782" ± 0.007" 
which is consistent with the astrometry of IDaemgen et all 
(2009) who found ARA = 0.763" ± 0.007", ADec = 
0.777" ± 0.007. The proper motion of TrES-2 is only 
HcCOsS — 2.34 zfc 1.7 mas/yr , ^t,s = —1.55 ± 1.7mas/yr 
(PPMX Catalog, iRoser et all 120081 ) . which is within our 
positional error bars over the time interval of ~ 2.4 years 
between observations. We have to await future observations 
to tell whether the pair is physically bound or not. 



3.2.3 TrES-4 

The candidate companion star was discovered 
in AstraLux Nort e obse rvations from June 2008 
by IDaemgen et al.l |2009). Our measurement of 
ARA = 0.003" ± 0.006", ADec = 1.550" ± 0.007" 
in October 2009 is consistent with the separation 
ARA = 005" ± 0.003" ADec = 1.555" ± 0.005" 
measured by IDaemgen efail (2009 ) . With a proper motion 
of Ha cos 5 = —9.9 4 ± 2.5 mas/yr, fis = —27.80 ± 2.5 mas/yr 
l|Roser et al.ll200sl ). the observations over a time baseline of 
16 months differ by ~ 1.4cr in RA and ~ 4.3<r in Dec from 
what would be expected if the companion candidate was a 
physically unrelated stationary background star, where a 
is the quadrature sum of the observed and proper motion 
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errors. We thus conclude that this is a common proper 
motion binary. 



companion, we can not exclude the possibility until future 
proper motion observations can be performed. 



3.2.4 WASP-2 

WASP-2 b was discovered bv lCollier Cameron et all l|2007l) . 
who also reported the stellar companion candidate to the 
east of the TEP host star at an angular separation of 
p = 0.7". The first AstraLux observation of this t arget 
was obtained in November 2007 (|Daemgen et alj|200gh . Al- 
though the primary star is a Kl V star, we list in Table [3] 
the derived magnitudes, colours and photometric distances 
assuming primary spectra l type K0 and K2, since co lours 
for Kl are not provided bv lKraus fc Hillenbrand! (|2007f ). De- 
pending on primary spectral type (see Sect. 3.1), we find that 
the companion candid ate is of spectral type M1/M2 ... M4. 
iDaemgen et~ail l|2009h observed the companion candidate at 
separation ARA = 0.732" ±0.004", ADec = 0.192" ±0.004" 
in November 2007. In November 2011 we observed the sepa- 
ration ARA = 0.720" ± 0.015", ADec = 0.188" ± 0.012". 
The proper motion of the planet host star is p a cos 5 — 
3.38 ± 2.9 mas/yr, [is = —52.31 ± 2.9mas/yr. The observa- 
tions over a time baseline of 4 years are marginally inconsis- 
tent with the hypothesis of a stationary background object 
in RA (~ 1.3<r) but by more than ller in Dec. We therefore 
conclude that this is a common proper motion system. 



3.2.5 WASP- 12 

The very b l oated planet WASP-12b was discovered by 
iHebb et al.l (|2009h . It is a highly irradiated planet, and 
one of the h ottest with an eq uilibrium temperature T cq = 
2516 ± 36 K ijHebb et al.ll2009l ). 

The stellar companion candidate is not previously 
known. We derive a spectral type K2 .. . M0 V from the mea- 
sured (i — z) colour and the SEDs of iKraus fc Hillenbrand! 
l|2007h . The proper motion of WASP- 12 is At a cos*5 = 
—0.3 6 ± 1.7 mas/yr, fig = —6.38 ± 1.7mas/yr (jRoser et alj 
2008), and future observations will be necessary to deter- 
mine whether or not the stars are physical companions. The 
companion candidate has an elongated shape (Fig. [2]) on 
both dates of observation, which suggests that the compan- 
ion candidate may be an unresolved binary. 



3.2.6 XO-3 

A very faint companion candidate (Az « 8.2) was found at 
large angular separation (p « 6") from the TEP host star 
XO-3. The (i — z)-colour places the companion candidate 
at a distance of > 860 pc if a main sequence star, and it 
is hence likely to be a non-related background object. An- 
other possibility is that the companion candidate is a white 
dwarf at the approximately same distance as the TEP host. 
The colours and brightness are consistent with a white dwarf 
with T e ji « 4000 K, corresponding to a very hydrogen defi- 
cient white d warf with cooling age 4 Gyr according to cool- 
ing curves bv lChen fc Hansen! i201lf). The age of X O-3 has 
been estimated to 2.821^2 Gvr (jWinn et al.ll2008h . and is 



3.3 Probability of chance alignment 

The probability of cha nce alignment is estim ated using the 
statistical approach of IDaemgen et ah! (|2009f ). The density 
of detectable background giants, p(m_R-), is c alculated by 
select ing all stars included in the 2MASS PSC (|Cutri et all 
12003( 1 within 30'of each of the observed targets (Tables [HE]) 
that are brighter than the estimated limiting magnitude of 
AstraLux, tuk ~ 14, and redder than (J — K) ^ 0.5. The 
probability of detecting a background giant is 



P(Q,m K ) 



-irp(m K )& 



(1) 



where 6 is th e maximum angular separation 
ijBrandner et al.l 120001 1. Using the aforementioned cuts 
in m K , (J - K) for 9 = 2" (see IDaemgen et~afll2009l ) the 
average probability of finding a non-related background 
star within 2"to the target star is P — 0.08%. We would 
then expect to detect unrelated background sources to 
E — 21 x P ~ 0.016 of our observed targets. This expec- 
tation value of chance alignment increases to 0.14 at the 
separation of 6" at which we find the probable background 
object in the observations of XO-3 (see Sect. 3.2.6). 

We also consider a possible contamination by M dwarfs 
by selecting stars in the 2MASS PSC with m K < 14.0 mag, 
and NIR colours J-H ^0.55 mag and H-K >0.0mag. We get 
a typical surface density of 6xl0 -5 M dwarfs per square 
arcsec, which corresponds to a probability of 7.5 xlO -4 to 
find an unrelated M dwarf within 2" of one of our target 
stars. We conclude that the probability of contamination by 
an unrelated M dwarf is of the same order as contamina- 
tion by a background giant. The close candidate compan- 
ions in our sample are thus likely to be true companions, 
although future observations are still necessary for confir- 
mation of common proper motions for HAT-P-8, TrES-2, 
HAT-P-7 and WASP-12. 



3.4 Bayesian statistical analysis 

We adap t the Bayesi an statistical analysis m ethod estab- 
lished in lAllenl l|2007h and utilized as well in IKraus et al.l 
(2011) to our survey for binarity of exoplanet host stars. 
We model the distributions of stellar binary mass ratios 
and projected separations as a power law in mass ratio 
(q = M scc /Mp r i m ), and a Gaussian in projected separation. 
Thus, we are left with 4 parameters to our models: 

• F - the companion star fraction (number of companions 
per primary star) 

• 7 - the power law index for the mass ratio (q) power 
law 

• ro - the average radius, i.e. the center of the Gaussian 
of projected separations 

• a - the width of the Gaussian of projected separations 



thus compatible with a 4.6 Gyr white dwarf at 3<r. While the 
candidate is unlikely to be a coeval comoving white dwarf 



The full combined form of our model is: 
R(q, r\F, 7, r a , a) oc Fq J exp(- 



: , (r-r a ) 2 



2a 2 



(2) 
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This equation gives us the probability of detecting a 
binary companion to a given star in a given projected sepa- 
ration, mass ratio range given a particular set of values for 
F, 7, ro and a. 

To calculate the likelihood, we must compare this model 
value with the actual data in each part of parameter space. 
To do so, we must calculate how many binary companions 
we expect to detect with this model in each projected sep- 
aration, mass ratio bin and then compare with the actual 
number detected in each bin (generally 0!). For a given mass 
ratio, projected separation bin and set of model parameters, 
the number of companions predicted will be: 

N pred (q, r) = N obs R(q, r\F, 7, r a , a) (3) 

where N j, s is the number of times this projected sepa- 
ration, mass ratio bin was observed in our survey (derived 
from the contrast curves and stellar properties of each sur- 
vey star). To compare data and model, we need to adopt a 
likelihood estimator. Since we expect to detect only small 
numbers of binary companions, our survey can be treated 
as a counting experiment. Thus, we adopt Poisson statistics 
here to calculate the likelihood: 

N^exp(-N det ) 
likelihood = prob{N det \F^,r ,o) = — — ■ (4) 

To derive the posterior probablitity distribution func- 
tion (PDF) for this bin, we must multiply the likelihood by 
any prior probability distribution for our parameters. Cer- 
tain priors can be adopted according to the functional form 
of the model parameters. We adopt the simplest uniform 
priors for 7 and r D : 

prob("/\I) — prob(r a \I) — 1 (5) 

However, F and a are scale parameters (i.e. invariant 
to changes in scale), allowing us to adopt a slightly more 
complex prior: 

prob(F\I) = i (6) 

prob(a\I) = - (7) 

a 

Thus our full prior is: 

prior — prob{F,^ ,r ,a\I) = — — (8) 

Fo 

Fo r a justification o f this choice, please see pages 109- 
110 in ISivia fc Skillind (120061 ). Multiplying the likelihood 
and prior then yields the posterior PDF for this projected 
separation, mass ratio bin. This can be generalized across 
all projected separation, mass ratio bins for the survey fairly 
easily. We generalize N j, s and Nrf e t into 2d arrays for each 
separation, mass ratio bin observed, which we will hence- 
forth call the window function and detection array. 

To build the window array, we use the contrast curve 
for each survey star to define the ranges in angular sepa- 
ration and contrasts where binary companions can be de- 
tected and the ranges where the contrast is insufficient to 
do so. Bins where a companion can be detected are assigned 
a value of 1; bins where no companion can be detected are 
assigned a value of 0. The window function for each star 
is then converted to projected separat ion vs. minimum de - 
tectable mass ratio using the models of lBaraffe et al.l (| 19981 ) 




-2-10 1 

y 

Figure 3. 2-D marginalized PDF for mass ratio power law index 
7 vs. companion star fraction F. 



and the distance to each survey star, z'-band contrasts and 
magnitudes were converted to H and K band for compar- 
ison with the iBaraffe et alj l| 19981 ) m odels by interpolating 
the st ellar SEDs from Table 5 from |K raus fc Hillenbrand! 
(2007). All the single star window functions are added to- 
gether to form the survey window function. The detection 
array is set up in a similar manner — first as a simple array 
with the number of objects detected in each separation and 
mass bin, and then deprojected into projected separation, 
mass ratio space. 

Then, for each set of model parameter, we calculate for 
each bin the posterior PDF. We multiply the posterior PDFs 
from each bin of observable space together to get the full 
posterior PDF across observable space for this set of model 
parameters. This process is repeated for all sets of model 
parameters of interest to derive the full posterior PDF as 
a function of the four model parameters. We calculated the 
posterior PDF for a grid in projected separation and mass 
ratio, running from 10 AU to 2000 AU in steps of 10 AU 
for projected separation and from 0.01 to 1.0 in q, in steps 
of 0.005. We allowed 7 to run from -3.9 to 3.7, in steps of 
0.4 and F to run from 0.05 to 9.72 in steps of 0.33. The 
average radius, r , ran from 10 AU to 2137 AU in steps of 
73.33 AU and a ran from 1 to 2901 AU in steps of 100 AU. 
Because it is not possible to visualize the full 4-dimensional 
posterior PDFs, we plot 2-d marginalized posterior PDFs in 
the 7— F and ro-a plane in Figures T3I4I The maximum of the 
posterior PDF yields the most-likely combination of model 
parameters given the information in hand. In this case, the 
maximum posterior value occurs with the parameter com- 
bination: r = 377 AU, o-=401 AU, F = 0.38, and 7=-0.7. 
However, from the marginalized posterior PDFs shown in 
Figures [3j [4] it is clear that many combinations of model 
parameters have fairly high likelihoods and that the covari- 
ance between parameters is fairly large. At the 67% (90%) 
confidence level, we constrain r to lie between 10 AU and 
890 AU (10 - 1476 AU), o to lie between 300 and 1800 AU 
(300 - 2901 AU), 7 to lie between -1.1 and 0.5 (-1.5 - 0.9) 
and F to lie between 0.38 and 2.05 (0.38 - 4.72). 
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Figure 4. 2-D marginalized PDF for the average projected bi- 
nary separation tq vs. the width of the Gaussian of projected 
separations a. 



rior PDF value of 7 = —0.7 (67% confidence), indicating 
a preference for low mass ratio and hence large magni- 
tude difference between companion stars. While this result 
is aga in in agreement with the field population of solar- type 
stars l|Duauennov fc Mavorll9 91; Ragh avan et al .l201fj|) .the 
blending by a nearby, bright star may cause rejection in tran- 
sit candidate follow ups, and thus the sample of TEP hosts 
may be biased towards faint, low mass-ratio companions. 

The main difference in multiplicity properties between 
field stars and planet host stars from this survey appears to 
be the average binary separation, which is about a factor of 
ten larger for the TEP hosts. With AstraLux we are typically 
sensitive to companions 2 magnitudes fainter than the planet 
host at angular separation of p = 0.2", corresponding to 
40 AU at the average TEP host distance of 200 pc in our 
target sample. All our discovered candidate companions are 
found at projected separations ^ 100 AU. While our sample 
is still small, the lack of close, faint companions may be a 
real property of exoplanet hosts. 



4 DISCUSSION 

The presence of a close binary companion is thought to affect 
processes of formation and subsequent dynamical evolution 
of the planets in such a system. We performed high angu- 
lar resolution Lucky Imaging of 21 host stars to transiting 
exoplanets in search of faint, close stellar companions. Two 
previously unknown companion candidates were discovered, 
to the TEP hosts HAT-P-8 and WASP-12. Future follow-up 
observations are necessary to confirm common proper mo- 
tion and hence physical companionship. Of the three candi- 
date binary TEP hos ts observed previously with AstraLux 
l|Daemgen et al.ll2009h . TrES-4 and WASP-2 were confirmed 
as common proper motion couples from follow-up observa- 
tions presented in this paper. Observations over a longer 
time baseline are still necessary to confirm the TrES-2 sys- 
tem as a physical binary, although we find that the prob- 
ability of non-related background objects is very low (see 
Sect. 3.3). A wide companion candidate to XO-3 is likely to 
be a background object, although we cannot rule out the 
possibility of a comoving white dwarf companion from these 
observations. 

We performed a Bayesian analysis of the full sample of 
31 TEP hosts o bserved with A s traLu x, including the sam- 
ple presented in iDaemgen et"aH (|2009h . The companion star 
fraction of the TEP host stars is n ot significantly differ- 
ent from that of solar type field stars (|Duquennov fe Mayor! 
ll99ll ; lRagh avan et al . 2010). This is in accordance with what 
was found in pr evious compilations of multiplicity among 
RV p l anet hosts (I Raghavan et ai1l2 006; Dcsi dera fe Barbieril 
120071 ; iBonavita fe Desideral 120071 ). At the 67% confidence 
level we find a companion star fraction lower limit of 
F ^ 0.38 whi ch is slightly h ig her t han the lower lim- 
its found by iRaghavan et al] (|2006l ) of w 26% and 
iMugrauer fe Neuhausei ( 20091 ). 20%. However, those mul- 
tiplicity fractions are based on samples of RV planet hosts 
and do not account for the common RV target selection bias 
of excluding known binaries from exoplanet surveys. No such 
pre-selection of single targets is made in transit surveys, and 
may explain the higher fraction derived in this survey. 

The slope of the mass-ratio distribution is likely to 
be uniform or slightly negative with a maximum poste- 
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